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ABSTRACT 

Identifying binaries among runaway O- and B-type stars offers valuable insight into the evo- 
lution of open clusters and close binary stars. Here we present a spectroscopic investigation of 12 
known or suspected binaries among field and runaway OB stars. We find new orbital solutions 
for five single-lined spectroscopic binaries (HD 1976, HD 14633, HD 15137, HD 37737, and HD 
52533), and we classify two stars thought to be binaries (HD 30614 and HD 188001) as single 
stars. In addition, we reinvestigate their runaway status using our new radial velocity data with 
the UCAC2 proper motion catalogs. Seven stars in our study appear to have been ejected from 
their birthplaces, and at least three of these runaways are spectroscopic binaries and are of great 
interest for future study. 

Subject headings: binaries: spectroscopic, stars: early-type, stars: kinematics, stars: individual (HD 
1976, HD 14633, HD 15137, Feige 25, HD 30614, HD 36576, HD 37737, HD 52266, HD 52533, HD 60848, 
HD 188001, HD 195592) 



1. Introduction 

Nearly all O- and B-type stars are believed to 
form in open clusters and stellar associations, but 
a handful of OB stars are observed at high galac- 
tic latitudes and with large peculiar space veloci- 
ties, suggesting that they have been ejected from 
the cluster of their birth. There are two accepted 
mechanisms to explain the origin of these runaway 
O- and B-type stars. In one scenario, close multi- 
body interactions in a dense cluster environment 
cause one or more stars to be scattered out of 
the region (Poveda, Ruiz, & Allen 1967). In rare 
cases, the ejected stars may be bound as a binary 
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pair; Leonard & Duncan (1990) use A^-body sim- 
ulations of open clusters to show that a binary 
frequency of about 10% is expected from dynami- 
cal ejections. An alternative mechanism involves a 
supernova explosion within a close binary, ejecting 
the secondary due to the conservation of momen- 
tum (Zwicky 1957; Blaauw 1961). The resulting 
neutron star may remain bound to the secondary 
if not enough mass is lost during the explosion. 
Portegies Zwart (2000) predicts a binary fraction 
of 20 — 40% among runaways which are ejected by 
such a binary supernova scenario. The observed 
fraction of binaries among runaways seems consis- 
tent with either scenario (5 — 26%; Mason et al. 
1998). 

Observational evidence suggests that both dy- 
namical interactions and supernovae in close bi- 
naries do produce some runaway systems. The 
microquasar LS 5039 was likely ejected from the 
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Milky Way disk about 1.1 Myr ago by the su- 
pernova explosion of the neutron star progenitor 
(Ribo et al. 2002; McSwain et al. 2004; Casares 
et al. 2005). Bergcr & Gics (2001) find that 
about 7% of Be stars have high peculiar space 
velocities, likely due to energy received during 
the supernova of a close binary companion (Mc- 
Swain & Gies 2005). On the other hand, Gua- 
landris, Portegies Zwart, & Eggleton (2004) show 
that a four-body encounter 2.5 Myr ago in the 
Trapezium cluster resulted in three ejected sys- 
tems: two single stars and one binary. The single 
stars both have very large space velocities and are 
classified as runaways, although the ejected binary 
has a lower space velocity due to its higher total 
mass. In a study of 44 trapczium-typc clusters, 
Allen, Poveda, & Hernandez- Alcantara (2006) find 
4 stars with large transverse velocities that were 
likely ejected by dynamical interactions. 

Which process dominates the production of 
runaway stars? Hoogerwerf, de Bruijnc, & de 
Zeeuw (2001) show that massive runaways tend to 
have enriched helium abundances and higher rota- 
tional velocities, while fewer non-runaways exhibit 
these traits. While rotational mixing could explain 
these results, it is doubtful that mixing occurs 
preferentially among runaways. It is more likely 
that prior mass transfer of CNO-processed mate- 
rial from the evolved donors has altered the sur- 
face abundances and added angular momentum to 
the mass gainers. The high binding energy of close 
massive binaries makes it unlikely they will be dis- 
rupted by dynamical interactions, but they could 
be separated by a supernova. Thus according to 
Hoogerwerf et al. (2001), the evidence favors the 
supernova ejection scenario over dynamical cluster 
interactions. However, pulsar searches by Philp et 
al. (1996) and Sayer, Nice, & Kaspi (1996) have 
not identified any neutron stars among runaway 
OB stars, and further investigation of their prop- 
erties is needed. 

The production of runaway OB binaries is ex- 
pected to be rare, but these systems can offer key 
insights into the evolution of close binary stars and 
open clusters. If massive stars form in dense re- 
gions through the mergers of low mass protostellar 
cores, the expected multi-body interactions will re- 
sult in a large number of high velocity runaways 
(Bally &: Zinnecker 2005). Ejecting a star dynam- 
ically from a young cluster will reduce the total 



energy imparted to the cluster via stellar winds, 
radiation, and supernova ejecta, whereas a proto- 
supernova system will remain in the cluster for a 
longer time and contribute more energy into its 
environment (Dray et al. 2005). Finally, the rela- 
tive importance of the supernova ejection scenario 
places limits on the mass needed to form black 
holes and neutron stars (Dray et al. 2005). 

In this work, we examine the binary status of 
12 field and runaway OB stars. Through a radial 
velocity study of these stars, we determine the or- 
bital solutions of five single-lined spectroscopic bi- 
naries (SBls). One is a possible triple star, and 
we also find two suspected SBl systems. By fitting 
the spectral energy distributions of these stars, we 
measure their reddening and angular sizes to con- 
strain their distances. Finally, our radial velocities 
combined with proper motions allow us to measure 
their peculiar spacx; velocities and reinvestigate the 
runaway status of these targets. Our results reveal 
that only 7 of the 12 stars have high space veloci- 
ties, while five stars may not be runaways. Of the 
runaways, two are SBls and a third is a suspected 
SBl. We will investigate their cluster ejection sce- 
narios in a forthcoming paper. 

2. Target Selection 

Most of our target stars (HD 14633, HD 15137, 
HD 30614, HD 37737, HD 52266, HD 52533, HD 

60848, HD 188001, and HD 195592) have been 
classified as field or runaway stars by Gies & 
Bolton (1986), Mason et al. (1998), and/or Stick- 
land & Lloyd (2001). Five of these (HD 30614, HD 
52266, HD 60848, HD 188001, and HD 195592) 
have unknown multiplicity or are suspected single- 
line spectroscopic binaries (SBls). The others are 
known SBl systems. Of this set, HD 14633 and 
HD 15137 are of particular interest since both SBl 
systems were likely ejected from the open cluster 
NGC 654, and their travel times appear longer 
than the expected O star lifetimes (Boyajian et 
al. 2005). On the other hand, the field stars HD 
52266, HD 52533, and HD 195592 may be the 
dominant members of previously undetected clus- 
ters (de Wit et al. 2004). 

HD 1976 has been classified as a member of the 
Cas-Tau OB association (De Zeeuw et al. 1999) 
and is a visual binary (Docobo & Costa 1986). It 
is also a slowly pulsating B star (Mathias et al. 
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2001) and a known SBl (Blaauw & van Albada 
1963). Although it is not a field star or a runaway, 
we include it in our observational program because 
its poorly known orbital solution (Blaauw & van 
Albada 1963; Abt et al. 1990) indicates a moder- 
ate eccentricity and a low mass function. Refining 
its orbit could help determine if this system is a 
post-supernova system whose kick velocity was too 
small to eject it from the association. 

The star HD 36576 is a field Bo star (Zorec, 
Fremat, & Cidale 2005) and a radial velocity vari- 
able (Gies & Bolton 1986). While Gies & Bolton 
(1986) classify this star as a non-runaway, Berger 
& Gies (2001) did not include this star in their 
sample when they performed a new investigation 
of high-velocity Be stars using Hipparcos proper 
motions. Therefore we reinvestigate its runaway 
status here. 

Feige 25 (HIP 12320) is a B-type star at high 

Galactic latitude (—48.36°) whose abundances 
and kinematics are consistent with a runaway 
Population I star (Martin 2004, 2006). Its ra- 
dial velocity has been measured only twice before 
by Martin (2006) and by Greenstein & Sargent 
(1974), and those measurements differ by nearly 
30 km s~^. Therefore we suspect that it may be 
an SBl system, and we include it here for further 
investigation. 

Most of the selected candidates are not known 
X-ray emitters, although three (HD 30614, HD 
52533, and HD 188001) are known weak X-ray 
sources. Mcurs, Fcnnell, & Norci (2005) showed 
that the observed X-ray luminosities and hardness 
ratios of HD 30614 and HD 188001 are consistent 
with the emission of normal O supergiant stars 
rather than from X-ray binaries. The hardness 
ratios of HD 52533 (Voges et al. 2000) are like- 
wise consistent with an origin in a normal O star 
(Motch et al. 1998). However, Meurs et al. (2005) 
point out that these results do not rule out the 
presence of a compact companion. If the accre- 
tion rate is sufficiently low, the resulting X-ray 
production may not contribute significantly to the 
total X-ray luminosity. A deeper X-ray observa- 
tion may be able to detect the characteristic hard 
spectrum of a MXRB in these systems. 

Likewise, none of our targets contain a known 
pulsar. Philp et al. (1996) searched for pulsars 
in 44 OB runaway stars, including our targets 
HD 30614, HD 37737, and HD 52533. Similarly, 



Sayer ct al. (1996) searched for pulsars in HD 
14633, HD 36576, and 38 other OB runaways. Al- 
though neither group found any pulsars associated 
with these stars, their existence cannot be ruled 
out. Using binary population synthesis calcula- 
tions, Portegies Zwart (2000) found that 20% - 
40% of OB runaways may have neutron star com- 
panions. However, due to the short lifetime of 
the pulsar (assumed to be 10 Myr) and the ab- 
sorption of the radio emission by the stellar winds 
(greater at periastron than apastron), Portegies 
Zwart (2000) predicts that only 1% - 2% of these 
neutron stars could be visible as radio pulsars at 
some time during their orbit. 

3. Observations and Radial Velocities 

Each target was observed during two observing 
campaigns at the KPNO 2.1 m telescope during 
2005 October and November. We used the #47 
grating (831 grooves mm~^) in 2"'' order to obtain 
a resolving power R = A/(5A ~ 3000. The observed 
wavelength range was 4050-4950 A, a region that 
includes numerous H-Balmer, He I, and He II lines 
in O- and B-type stars. The data were reduced 
using standard routines in IRAF. 

We also obtained some spectra from the KPNO 
0.9 m coude feed telescope in 2005 November. The 
resolving power is R ~ 12000 using the long colli- 
mator, grating A (in second order with order sort- 
ing filter 4-96), camera 5, and the T2KB CCD, a 
2048 X 2048 device. We obtained a spectral cov- 
erage of 4240-4580 A. 

We used all available strong, unblended lines 
to measure radial velocities, Vr, in these stars. 
Rest wavelengths for each line were taken from the 
NIST Atomic Spectra Databasc-'^. We selected one 
spectrum with good signal-to- noise (S/N) to serve 
as a template, and we fit the core of each absorp- 
tion line with a parabola to determine its absolute 
radial velocity. The remaining spectra were then 
cross correlated with this template to determine 
the relative velocities. Several lines in the 2.1 m 
data sets lie near bad pixels which corrupted the 
cross correlations, so these lines were usually omit- 
ted. In a few cases we found systematic difFc^rc^nces 
between sets of lines, and we discuss these below. 
Only in the case of HD 1976, we found a system- 

The NIST Atomic Spectra Database is available online at 
http://physics.nist.gov/PhysRefData/ASD/index.html. 
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atic disagreement between our measurements and 
those in the hterature. Our Vr measurements of 
all other stars were in good agreement with pre- 
vious measurements. For each star, we list the 
complete set of lines used in Table 1, and the re- 
sulting mean velocities are listed in Tables 4 and 
5. For the stars in Table 5, we also include the 
standard deviation, a, of the Vr measured from 
many lines for each observation. 

To test the significance of systematic night-to- 
night (NN) and line-to-line (LL) variations, we 
performed a 2-way analysis of variance (2A0V) 
statistical test for each set of velocities. This test 
is described in detail by Gies & Bolton (1986). 
Essentially, it allows us to test the null hypothe- 
sis that the star does not exhibit statistically sig- 
nificant variations. The F statistic describes the 
variance of the entire data set divided by the av- 
erage variance of each measurement, and p gives 
the probability that the observed variations are 
drawn from the same random distribution. If p 
is small, less than 1%, then we must reject the 
null hypothesis. Thus a low p indicates that the 
observed differences in Vr are significant. Both F 
and p depend upon the number of degrees of free- 
dom, which is just the number of input variables 
(lines or spectra) minus one. 

The coude feed spectral range is limited and in- 
cludes only a few of the lines measured in the 2.1m 
spectra. Consequently, we applied the 2A0V test 
to the common line sample in 2.1 m data alone, 
and these results are summarized in Table 2. The 
test indicates significant LL variations in HD 1976, 
HD 15137, HD 30614, HD 37737, HD 52533, and 
HD 60848. We are not surprised to find LL vari- 
ations in OB-type stellar spectra because line for- 
mation may occur at different levels in an expand- 
ing atmosphere and because some lines sniffer from 
blending with weaker nearby components. Fur- 
thermore, HD 52533 may have an optical com- 
panion causing contamination in its line profiles. 
Even though these stars do have statistically sig- 
nificant LL variations, we can still identify signif- 
icant nightly variations in Vr since we measure a 
consistent set of lines in the our spectra. Feige 
25 and HD 60848 show little evidence of nightly 
systematic differences, while the nightly variations 
observed in all other targets are highly significant. 

For the stars with significant night-to-night Vr 
variations, we combined data from the literature 



(where available) with our own measurements for 
improved orbital solutions. We performed period 
searches using a version of the discrete Fourier 
transform and CLEAN dcconvolution algorithm 
of Roberts, Lehar, & Dreher (1987) (written in 
IDL^ by A. W. Fullerton). We then tested any 
significant signal found in the CLEANed power 
spectrum to verify that it produced a reasonable 
phased Vr curve. We ruled out periods less than 
one day because we never observed rapid, high- 
amplitude changes in Vr among our targets. The 
best period from the CLEANed power spectrum 
was used with the non-linear, least-squares fitting 
program of Morbey & Brosterhus (1974) to solve 
for the orbital elements. 

A detailed discussion of the results from our ra- 
dial velocity study is presented below, and we sum- 
marize the results here. We were able to improve 
the known orbital elements for 5 SBl systems (HD 
1976, HD 14633, HD 15137, HD 37737, and HD 
52533), and these are listed in Table 3. Their Vr 
curves are shown in Figure 1 , and our Vr measure- 
ments are listed in Table 4. The deviations from 
the theoretical Vr curves are comparable to the 
errors of the mean. 

The stars Feige 25, HD 36576, and HD 60848 
appear to be constant velocity stars. Both HD 
36576 and HD 60848 are Be-type objects, so their 
Vr are affected by variable emission filling the 
line profiles or photospheric nonradial pulsations 
changing their line shapes. HD 30614 and HD 
188001 have been identified as SBl systems in the 
past, but we use available data in the literature 
to show that they are more likely single stars that 
have been misclassificd. HD 52266 and HD 195592 
arc likely SBl systems, but we cannot determine 
the systems' periods or orbital elements from the 
available data. The Vr measurements for these 
stars without orbital solutions in this work are 
listed in Table 5. 

3.1. HD 1976 

HD 1976 is a speckle binary with SV = 0.89 
and separation 0'.'138 (Hartkopf et al. 2000; Ma- 
son et al. 2001). It is also an SBl, and Blaauw 
k van Albada (1963) and Abt et al. (1990) each 
present orbital solutions for HD 1976. These solu- 
tions have somewhat differing orbital period (27.8 

^IDL is a registered trademark of Research Systems, Inc. 
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d and 25.44 d, respectively) and eccentricity (0.2 
and 0.14, respectively). 

Unfortunately, our observing schedule did not 
allow us to obtain good orbital coverage of this sys- 
tem, and we rely heavily upon the data of Blaauw 
& van Albada (1963) and Abt et al. (1990) to de- 
termine the period and orbital elements. To up- 
date the orbital elements for this SBl, we included 
our own 15 Vr measurements with those published 
by Blaauw & van Albada (1963) and Abt ct al. 
(1990) (except 4 points with large scatter, which 
we assigned half weight). We noticed a system- 
atic difference between their data and ours, so we 
added an offset of 4.5 km to the velocities of 
Blaauw & van Albada (1963) and 5.3 km s~^ to 
the velocities of Abt et al. (1990) to bring our 
measurements into better agreement. Our period 
search revealed two possible results, 25.4 or 27.6 
d, and we investigated many possible periods close 
to these values to obtain the lowest errors in the 
orbital fit. The resulting errors were two times 
smaller using the shorter period, and we adopt a 
period of 25.4176 d. The remaining orbital ele- 
ments for HD 1976 are listed in Table 3. 

We expect that the speckle companion con- 
tributes about 30% of the flux in this system, but 
we saw no signs of blending in the H Balmer or 
He I line profiles. However, the relatively low ve- 
locity semiamplitude, K, and high V smi of this 
star (see Section 4) may prevent us from detecting 
the spectrum of the companion star. 

3.2. HD 14633 

In our earlier Vr study of HD 14633 (Boyajian 
et al. 2005), we determined that this SBl has a 
short period (15.4083 ± 0.0004 d), high eccentric- 
ity, and low mass function. This system may har- 
bor a compact companion formed by a past su- 
pernova, making it one of the first known quiet 
MXRB candidates. In Table 4, we present 14 new 
Vr measurements of this system. 

To improve the orbital period, we combined 
a total of 103 Vr measurements from Bolton 
& Rogers (1978), Stone (1982), Boyajian et al. 
(2005), and this work. We refine the period to 
15.407 d. From the recent data of Boyajian et al. 
(2005) and this work, we also improved the orbital 
elements, listed in Table 3. One point at cj) ~ 0.4 
has a very large scatter, and we assigned it zero 



weight in the fit. The new orbital elements are the 
same within error as our earlier values (Boyajian 
et al. 2005). 

3.3. HD 15137 

Likewise, HD 15137 may also be a quiet MXRB 
candidate (Boyajian et al. 2005) which we ob- 
served to refine its orbital elements. We measured 
Vr for HD 15137 using the same Balmer, He I, and 
He II lines we used for HD 14633. Because of the 
broad lines in the spectrum of HD 15137, we find 
significantly more scatter than in the case of HD 
14633. Our Vr measurements are listed in Table 
4. 

For HD 15137, we collected 34 published V. 
measurements from Conti, Lecp, & Lorre (1977) 
and Boyajian et al. (2005) to include with our 14 
new measurements. A period search using this 
complete data set suggests a period of 30.35 d, 
slightly higher than our preliminary result (28.61 
d; Boyajian et al. 2005). The addition of our 
new measurements also affects the other orbital 
elements, and our revised orbital solution is pre- 
sented in Table 3. However, we stress that the 
period of HD 15137 relies heavily on the single 
observation from Conti et al. (1977) and thus it 
remains highly uncertain. Better orbital coverage 
will improve the results significantly. 

3.4. Feige 25 

On five of the seven nights we observed Feige 
25, we obtained two 30 minute exposures which 
wc coadded for better S/N. On the remaining 
nights, we were only able to obtain a single 30 
minute exposure. The high a in our measure- 
ments likely results from the low S /N in our spec- 
tra since the noisy He I and Mg II lines were diffi- 
cult to fit. While we observe some scatter in the Vr 
set, it is not significant and the velocity remains 
constant within the measured errors. Since our 
measurements are also consistent with the values 
Vr = 24.1 ± 12.1 (Martin 2006) and Vr = -5 ± 25 
(Greenstein & Sargent 1974), we classify this star 
as radial velocity constant. 

3.5. HD 30614 

HD 30614 has been classified as an SBl by Zej- 
nalov & Musaev (1986), who found a 3.68 d period 
and highly eccentric (e = 0.45) orbit with a very 
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low velocity scmiamplitudc, K ~ 9 km s~^. Such a 
low K would imply a very low inclination in a mas- 
sive SBl, and selection effects make it extremely 
difficult to detect such systems. We combined our 
15 observations with additional measurements of 
its Vr from Beardsley (1969), Conti et al. (1977), 
Bohannan & Garmany (1978), Stone (1982), Gies 
& Bolton (1986), and Zejnalov & Musaev (1986). 
A new period search using this complete data set 
revealed the possible P = 3.24 d. However, these 
studies each used the mean Vr from many lines, 
and the interline scatter is comparable to the low 
K fomid by Zejnalov & Musaev (1986). 

Instead of using all of these data, we chose to 
investigate a small subset of lines with the hope of 
reducing the scatter in the radial velocity curve. 
Only Zejnalov & Musaev (1986) and Gies (private 
communication) provide Vr measurements for in- 
dividual lines, and three lines are common to our 
data sets: Si IV AA4089, 4116 and B.5. Since HD 
30614 is a non-radial pulsator with variable stel- 
lar winds (FuUerton, Gies, & Bolton 1996; Rzayev 
& Panchuk 2004), the Balmer lines may be con- 
taminated. Therefore we chose to use the mean 
from only the two Si IV lines, which usually agree 
within 4 km s~^ in our data. We repeated the pe- 
riod search using this subset of data, and we found 
a strong signal at P 3.57 d, closer to that found 
by Zejnalov & Musaev (1986), but the low ampli- 
tude variations probably originate in atmospheric 
or wind modulations rather than orbital motion. 
Without the availability of higher resolution spec- 
tra, we conclude that HD 30614 is more likely a 
single star that has been misclassified. 

3.6. HD 36576 

Because the H Balmer lines have strong emis- 
sion in this B2 IV- Ve star, we used only the avail- 
able He I lines to measure Vr, listed in Table 5. 
Even among these lines, their asymmetric profiles 
caused some difficulty in our measurements, re- 
sulting in a large scatter for each Vr- The 2A0V 
test confirms that the high a is due to our mea- 
surement errors rather than real velocity differ- 
ences between lines. In our sequence of observa- 
tions, we find rapid yet low-amplitude variations 
in Vr that suggest non-radial pulsations may be 
present in this star. 



3.7. HD 37737 

An orbital solution for HD 37737 was originally 
published by Gies & Bolton (1986) with an orbital 
period of 2.49 d. However, our period search using 
our data with theirs resulted in a very different 
period of 7.84 d. Even using only Gies & Bolton's 
data, we found the same period with the modern 
CLEAN algorithm. We conclude that the original 
period was in error, and we present new orbital 
elements for this SBl in Table 3. Figure 1 shows 
that the Gies & Bolton (1986) data agree well with 
our data and the new orbital solution. 

3.8. HD 52266 

We measured Vr for HD 52266, listed in Ta- 
ble 5, using the mean from many H-Balmer, He I, 
and He II lines. The star has very broad lines 
that contribute to the high scatter in our mea- 
surements, and the short term changes in Vr are 
usually smaller than the measured errors. How- 
ever, we find a definite increase in Vr between the 
October and November runs that is statistically 
significant in our 2A0V test. We suspect that HD 
52266 may be a SBl with a period too long to 
measure from our data. 

3.9. HD 52533 

We initially measured Vr for HD 52533 using 
many H Balmer, He I, and He II lines. The mean 
velocity from the Balmer lines varies from 3 — 122 
km s""'^, and comparable variations are seen in the 
He I lines. However, the He II line measurements 
have a much higher amplitude of variation, from 
— 18 to 206 km s"-'^! These lines indicate significant 
line-to-line variations in this star, and the 2A0V 
test reveals that the probability of these varia- 
tions originating from random fiuctuations within 
the same sample is negligible. We suspect that 
this may be a triple star system; the He II line 
probably originates in the hottest component, an 
0-type primary star in a close spectroscopic bi- 
nary. A cooler, B-type companion is probably a 
widely separated, stationary star. Its strong He I 
and Balmer lines cause blending in these line pro- 
files without affecting the He II lines. HD 52533 
is not a known speckle binary. Most recently, it 
was observed but not resolved by B. Mason in 
2005 November (private communication), placing 
an upper limit of 30 mas on the angular separation 
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of the wide components assuming SV < 3. 

Because of the suspected Une blending, we in- 
vestigate the orbit of the primary using the mean 
Vr from only the He II AA4200, 4542, 4686 lines. 
Gies (private communication) provided 10 mea- 
surements of these lines that we combine with our 
20 observations to determine the primary's orbital 
period, P = 22.19 d. Only one He II line was 
recorded in our coude feed spectra, and the large 
scatter in Gies's data led us to use only one or two 
of his He II line measurements in some cases. For 
each Vr based on only one line instead of a mean, 
we assigned the point only half weight. The re- 
sulting orbit of the primary is presented in Table 
3, and its Vr measurements are listed in Table 4. 
We will explore this system in more detail in a 
forthcoming paper. 

3.10. HD 60848 

The Bf], H7, and He I A4713 Uncs in HD 60848 
show double peaked emission in our spectra, there- 
fore we omitted these lines in our Vr measure- 
ments. Like HD 36576, we found rapid, low am- 
plitude variations in Vr in HD 60848 that are con- 
sistent with non-radial pulsations, and we classify 
this star as single. The high scatter in our mea- 
surements is not statistically significant, and it is 
likely due to the pulsations' effect on line profiles 
in HD 60848. 

3.11. HD 188001 

There are a wealth of Vr measurements for 
HD 188001 in the literature, and several groups 
(Underbill & Matthews 1995; Aslanov, Kornilova, 
& Cherepashchuk 1984; Aslanov & Barannikov 
1992) have published orbital solutions for this 
spectroscopic binary with periods ranging from 
32.5—78.74 d. Our new observations do not agree 
with any of the previous solutions, so we decided 
to undertake a new analysis of this star using the 
available data. Thus we gathered 14 measure- 
ments of Vr from Plaskctt (1924), 2 from Conti et 
al. (1977), 11 from Bohannan & Garmany (1978), 
8 from Garmany, Conti, & Massey (1980), 16 from 
Stone (1982), 17 from Aslanov ct al. (1984), 20 
from Fullerton (1990), 9 from Underbill (1995), 
and 8 new values from this work (listed in Table 
5). 

This complete sample of 105 values is an in- 



homogeneous collection of measurements relying 
on different line groups and different wavelength 
regions, but the measured lines are predominantly 
He I and He II with some metal lines also included. 
To retain the maximum consistency between these 
data sets, whenever the authors present Vr from 
more than one line or line group, we chose to use 
only He I and He II lines in absorption. These lines 
also provided the most reliable Vr in our own mea- 
surements; however, they are probably affected by 
the stellar winds. We observe the He II A4686 line 
in emission, and Fullerton (1990) notes a variable 
P Cygni absorption trough in the He I A5876 line 
profiles, whereas the C IV AA5801, 5812 lines in 
Fullerton's sample are more likely photospheric. 
Certainly we arc concerned by the need to rely on 
wind lines to investigate the orbital period of HD 
188001, but in lieu of a more complete data set 
from purely photospheric lines, we are limited in 
our approach. 

We performed a period search with these 105 
values, eliminating spurious frequencies that sug- 
gest periods less than 10 days since rapid Vr 
variations are not observed in HD 188001. The 
strongest signal suggests a period P = 29.83 d; 
however, the sinusoidal variation appears to be 
due to the systematic differences between various 
data sets rather than true orbital motion. There- 
fore we conclude that HD 188001 is actually a sin- 
gle star, possibly with variable winds that con- 
taminate the He absorption lines and mimic the 
signature of a spectroscopic binary. 

3.12. HD 195592 

Wc did not use the H Balmcr lines to measure 
Vr in HD 195592 due to their inverse P Cygni 
profiles, suggesting that wind emission is partially 
filling in these lines. Instead wc used strong, un- 
blended He I, He II, and Si IV lines to measure 
velocities. We find that Vr is slightly lower dur- 
ing our November nm than in October, and the 
2A0V statistical test indicates that this differ- 
ence is highly significant. Therefore we believe 
HD 195592 may be a SBl system. Wc combined 
our data with 8 Vr measurements from Mayer et 
al. (1994) to search for a possible orbital period, 
but the only significant signal occurcd at the un- 
physical period ~ 0.5 d. Additional data is needed 
to determine if HD 195592 is indeed a binary. 
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4. Spectroscopic Modeling and SED Fits 

Our O- and B-type targets include a wide range 
of spectral types, luminosity types, and even emis- 
sion stars, making it difficult to obtain the temper- 
ature and surface gravity of every star. Nonethe- 
less, we used a combination of sources to obtain 
their physical parameters. The resulting values, 
and references where applicable, are listed in Ta- 
bic 6. 

For the B-type stars HD 1976 and Feige 25, we 
generated a grid of synthetic, plane-parallel, local 
thermodynamic equilibrium (LTE) atmospheric 
models using the Kurucz ATLAS9 code (Kurucz 
1994). We adopted solar abundances and a mi- 
croturbulcnt velocity of 2 km for these stars, 
which corresponds to the mean microturbulence 
observed among late-type, main-sequence B stars 
(Lyubimkov. Rostopchin, & Lambert 2004). Each 
atmospheric model was then used to calculate a 
grid of model spectra using SYNSPEC (Lanz & 
Hubeny 2003). Using the measured Vr from each 
observation, we shifted each observed spectrum to 
its rest frame and created a mean spectrum of each 
target to compare with the model grid. We com- 
pared the observed three H Balmer line profiles 
to model profiles convolved with a limb- darkened, 
rotational broadening function and a Gaussian in- 
strumental broadening function to measure T^g, 
logg, and Vsini. For each fit, we used steps of 
100 K in Toff, 0.1 dex in logg, and 10 km in 
Vsini and determined the best fit by minimizing 
the square of the residual, (O — C)^, over a 40 
A wide region centered on the rest wavelength of 
each line. The line wings are particularly sensitive 
to loggf while the line depths are strong indica- 
tors of Toff since the strength of the Balmer lines 
declines steadily with increasing T^g in the B spec- 
tral regime. Finally, the shape of the line core is 
most sensitive to Vsini, so we were able to obtain 
excellent fits by adjusting these three parameters. 
Using this technique, we find errors of ±100 K 
in Teff for Feige 25 and ±400 K for HD 1976. For 
both stars, the error in logg is ±0.1 and for V sini 
it is ±30 km s^^. The resulting parameters for HD 
1976 and Feige 25 are listed in columns 3 — 5 of Ta- 
ble 6. Our measurements of Feige 25 are in good 
agreement with the physical parameters found by 
Keenan & Dufton (1983) and Martin (2004). 

For the hotter 0-type stars, we used the OS- 



TAR2002 grid of line-blanketed, non-LTE, plane- 
parallel, hydrostatic atmosphere model spectra 
for 0-type stars from the TLUSTY code (Lanz & 
Hubeny 2003). These models adopt solar abun- 
dances and a fixed microturbulent velocity of 10 
km s~^. Once again, we convolved each model 
with a limb darkened rotational broadening func- 
tion and a Gaussian instrumental broadening 
function to compare with our observed spectra. 
Here we used step sizes of 100 K in Toff-, 0.05 in 
logg, and 10 km in Vsini. Because of the 
probable line blending in the Balmer and He I 
lines in HD 52533 and emission partially filling 
the He I line profiles in HD 60848, we used the 
He 11 AA4200, 4542, 4686 lines to measure the 
physical parameters of each 0-typc dwarf. The 
resulting T^g, logg, and Vsini from the O star 
fits are listed in Table 6. Our measurement errors 
are typically 400 K in T^g, 0.10 in logg, and 10 
km s~^ in Vsini. 

The Be star HD 36576 had obvious emission in 
its Balmer lines, and the O- supergiants HD 30614, 
HD 188001, and HD 195592 were also difficult to 
fit with the plane-parallel spectral models. There- 
fore we adopted Toff and log g of HD 30614 and HD 
36576 from Crowther, Lennon, & Walborn (2006) 
and Chauville et al. (2001), respectively. In ad- 
dition, we adopted Vsini = 165 km for HD 
36576 (Abt, Levato, k Grosso 2002). We used 
the spectral types of HD 188001 and HD 195592 
to obtain their parameters from the calibration of 
Martins, Schaerer, & HiUier (2005). 

We then compared the resulting model spec- 
trum to the observed spectral energy distribution 
(SED). For all of our targets except Feige 25 and 
HD 52533, flux-calibratcd UV spectra were avail- 
able from the International Ultraviolet Explorer 
(lUE) archives. We also obtained UV fluxes in 
several bandpasses from Thompson et al. (1978) 
and/or Wesselius et al. (1982) in most cases. For 
HD 52266 and HD 195592, the available WE 
fluxes disagree with the other UV fluxes, so we pre- 
fer the UV bandpass photometry which provides 
a continuous SED when combined with optical 
data. Johnson U BVR (from numerous sources in 
the VizieR database; Ochsenbein, Bauer, & Mar- 
cout 2000), Stromgren uvby (Hauck & Mermilliod 
1998), and 2MASS JHK (Cutri et al. 2003) pho- 
tometry were usually available as well. For these 
three fllter systems, we converted the photometric 
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magnitudes to fluxes using the techniques of Col- 
ina, Bohlin, & Castelli (1996), Gray (1998), and 
Cohen, Wheaton, fc Megeath (2003), respectively. 
The model spectra were binned into 50 A bins in 
order to remove small scale line structure and com- 
pare to the observed SED. For those stars with 
close binary or speckle companions that contribute 
to the SED, we neglect the companion's flux con- 
tribution in the model. 

Using a grid of values for the reddening, E{B — 
V), and the ratio of total-to-selective extinction, 
i?, we compared the reddened absolute fluxes from 
the spectral models to the observed stellar flux. 
We used a step size of 0.01 for both parameters, 
and we determined the best fit values of E{B — V) 
and R by minimizing (O — C)^, and the typical er- 
rors are 0.03 mag and 0.4, respectively. The ratio 
of the observed stellar flux to the reddened model 
flux provides the angular diameter, ^lDj of each 
star (Gray 1992). The errors in 9ld are deter- 
mined from the standard deviation of the ratios 
from each optical and infrared flux estimate since 
the ratio is wavelength independent. These mea- 
sured parameters are listed in columns 6 — 8 of 
Table 6. 

Distances to most of our targets were obtained 
from van Steenberg & ShuU (1988), Mason et al. 
(1998), Patriarchi, MorbidelU, & Perinotto (2003), 
and Martin (2006). For the remaining stars HD 
1976 and HD 36576, we obtained an estimate 
of their absolute magnitudes from Wegncr (2000) 
and combined these with the observed V and our 
calculated A{V) = R x E{B — V) to determine 
their distances. Finally, we used Old with the 
known distances to determine the radius, R^,, for 
each star (column 9 of Table 6). We note that 
the closest stars in our list, HD 1976 and HD 
36576, also have Hipparcos parallaxes available, 
but the larger distances from these measurements 
(420 and 570 pc, respectively) imply significantly 
larger R^, than we expect for these stars. 

In Table 6 we also provide the predicted i?* 
from the calibration of Martins et al. (2005) for 
stars with T^if > 30000, and for cooler stars we use 
the evolutionary models of Schaller et al. (1992) to 
determine the expected value. The observed dif- 
ferences between the derived and predicted radii 
are usually small, and we find no systematic dif- 
ferences between the binaries' and single stars' 
radii that may be due to neglecting the compan- 



ion fluxes in the SED fits. Note that wo do find 
large discrepancies in R^, for the binaries HD 37737 
and HD 195592, but we find a similar error for 
the single star HD 188001. All three stars are gi- 
ants and supergiants, and the errors in their pre- 
dicted physical parameters are larger (Martins et 
al. 2005). In these cases the distances require a sig- 
nificant downward revision to match our derived 
radii with the expected values. Therefore we de- 
rive a new distance to each target based upon the 
expected Ri, and our measured 0ld , and we adopt 
the mean of these two distances with errors rep- 
resenting the range of possible values. The final 
adopted distances are listed in Table 6. 

We also inspected the observed SED to look for 
an IR excess in our targets, which may be indica- 
tive of an accretion disk around a compact com- 
panion. The emission-line stars HD 36576 and HD 
60848 appear to have a small iiT-band flux excess 
due to their circumstellar disks, but no other tar- 
gets show evidence of an IR excess. 

5. Runaway Status of Targets 

Our new radial velocity data and the distance 
constraints above allow us to reexamine the run- 
away status of each of our targets. We list in Table 
7, columns 2 — 3 the galactic longitude. £, and lat- 
itude, b, for each star. Columns 4 — 5 give the 
proper motions, /x„cos(5 and yU^, obtained from 
The Second U.S. Naval Observatory CCD Astro- 
graph Catalog (UCAC2; Zacharias et al. 2004) 
and The UCAC2 Bright Star Supplement (Urban, 
Zacharias, & Wycoff 2004). We computed the pe- 
culiar transverse and radial space velocity, Vt^pec 
and Vr.pec , of each star using the method of Berger 
& Gies (2001), and these are listed in columns 
6 — 7. These are combined into a single measure 
of the peculiar space velocity, Vpec, in column 8. 
Our values of Vt^pec and Vr^pec for HD 37737 and 
HD 52533 agree closely with those found by Moffat 
et al. (1998). 

We classify as definite runaways those stars 
with Vpec > 30 km s~^ plus the error in Vpec 
(Gies & Bolton 1986), and we find seven such stars 
among our sample: HD 14633, HD 15137, Feige 
25, HD 30614, HD 36576, HD 188001, and HD 
195592. Our findings agree with those of Noriega- 
Crespo, van Buren, & Dgani (1997), who found 
bow shocks associated with HD 30614, HD 188001, 
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and HD 195592 due to their fast space velocities. 
Uncertainties in the proper motions and distances 
contribute to large errors in Vpec for several of our 
stars, so the tentative runaways HD 37737 and HD 
52533 also have Vpec > 30 km s~^ within their er- 
rors. Having a lower Vpec does not exclude the 
possiblity that a star was ejected in the past as 
long as its velocity is greater than the escape veloc- 
ity from the cluster, an ejected star will effectively 
be a runaway. 

Although we find that HD 37737 is tentatively a 
runaway, it is worthwhile to note that the star has 
also been classified as a member of the Aur OBI 
association (Humphreys 1978). It may also be 
surrounded by a symmetric H II bubble (Noriega- 
Crespo et al. 1997), although the morphology of 
the bubble is difficult to classify with certainty. 
To investigate whether HD 37737 is a member 
of Aur OBI, we compared its proper motion to 
those of the other proposed Aur OBI members 
using the UCAC2 and UCAC2 Bright Star Sup- 
plement catalogs. We find that HD 37737 and 7 
other stars from Humphrey's list do indeed have 
similar proper motions, although the remaining 7 
stars are moving in nearly an opposite direction. 
From these results we doubt that the Aur OBI as- 
sociation is truely a bound group, but the evidence 
does suggest that multiple open clusters may be 
superimposed upon this region of the sky. Further 
investigation of the groups' distances are required 
to determine whether HD 37737 is a member. 

The stars HD 1976, HD 52266, and HD 60848 
have low Vpec and are probably not runaway stars. 
HD 52266 (along with the higher velocity systems 
HD 52533 and HD 195592; de Wit et al. 2004) may 
be a member of a previously undetected cluster. 
Since our spectra of HD 52533 suggest that it may 
be a triple system, it is very likely that a sparse 
cluster indeed exists and that this system is not 
a runaway. We will explore this possibility in the 
future. 

Our search was designed to detect runaway bi- 
nary systems, and it is no surprise that we iden- 
tify several in our sample (HD 14633, HD 15137, 
HD 195592, and the possible runaway HD 37737). 
Both dynamical interactions within a dense open 
cluster and a supernova in a close binary can pro- 
duce runaway SBl systems, and our results high- 
light how rare such systems are. These ejection 
mechanisms can be distinguished observationally, 



and we explore both scenarios for HD 14633 and 
HD 15137 in a forthcoming paper. 
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Table 1 

Lines Used for Radial Velocities 



Star 


Telescope 


Lines 


HD 1976 


2.1 m 
coude 


H/3, H7, H(5, He I AA4144, 4388, 4471, 4713, 4922, C II A4267, Mg II A4481 
H7, He I AA4388, 4471 


HD 14633 


2.1 m 


H/3, H7, US; He I AA4121, 4144, 4388, 4471, 4713 (Oct. only), 4922; He II AA4200, 4542, 4686; Si IV A4116; N III A4640* 


HD 15137 


2.1 m 


H/3, H7, HS, He I AA4144, 4388, 4471, 4713, 4922; He II AA4200, 4542, 4686; Si IV A4089; N III A4379 


Feige 25 


2.1 m 
coude 


H/3, H7, HS; He I AA4144, 4388, 4471, 4922 (Oct. only); Mg II A4481 
H7; He I AA4388, 4471 


HD 30614 


2.1 m 

coude 


H/3, H7, HS; He I AA4121, 4144, 4388, 4471, 4713, 4921; He II AA4200, 4542; Si IV AA4089, 4116; N III A4379 
HS, Ho I AA4388, 4471; Ho II A4542; N III A4379 


HD 36576 


2.1 m 
coude 


He I AA4121, 4144, 4388, 4471, 4713 
He I AA4388, 4471 


HD 37737 


2.1 m 

coude 


H/3, i?7, HS; He I AA4144, 4388, 4471, 4713, 4922; He II AA4200, 4542, 4686; Si IV A4089 
H-y; He I AA4388, 4471; He II A 4542 


HD 52266 


2.1 m 

coude 


H/3, Hj, HS; He I AA4144, 4388, 4471, 4713, 4921; He II AA4200, 4542, 4686; Si IV A4089 

i?7; Ho I AA4388, 4471 


HD 52533 


2.1 m 
coude 


He II AA4200, 4542, 4686 
He II A4542 


HD 60848 


2.1 m 
coude 


HS; Ho I AA4144, 4388, 4471, 4921; He II AA4200, 4542, 4686; Si IV A4089; N III A4379 
He I AA4388, 4471; He II A4542 


HD 188001 


2.1 m 

coude 


He I AA 4120, 4144, 4388, 4471, 4713, 4922; He II AA4200, 4542 
He I AA4388, 4471; He II A4542 


HD 195592 


2.1 m 

coude 


He I AA4120, 4143, 4387, 4471, 4713, 4922; He II AA4200, 4542; Si IV AA4089, 4116; N III AA4379, 4640* 
He I AA4388, 4471; He II A 4542 



"^Blend. 



Table 2 

Results from Two-way Analysis of Variance Test 



Star 


Degrees of 
Freedom (LL) 


F(LL) 


P(LL) (%) 


Degrees of 
Freedom (NN) 


F(NN) 


p(NN) (%) 


HD 1976 


10 


2.54 


0.78 


13 


13.45 


0.00 


HD 14633 


13 


1.96 


2.72 


13 


33.65 


0.00 


HD 15137 


12 


5.17 


0.00 


13 


3.18 


0.03 


Feige 25 


7 


1.43 


23.32 


4 


2.74 


4.84 


HD 30614 


13 


11.07 


0.00 


13 


12.17 


0.00 


HD 36576 


4 


2.08 


9.87 


12 


8.77 


0.00 


HD 37737 


11 


14.39 


0.00 


12 


316.90 


0.00 


HD 52266 


11 


1.84 


5.22 


14 


10.38 


0.00 


HD 52533 


9 


6.48 


0.00 


14 


13.30 


0.00 


HD 60848 


9 


11.35 


0.00 


13 


2.22 


1.27 


HD 188001 


7 


1.54 


18.09 


6 


11.36 


0.00 


HD 195592 


11 


2.05 


3.41 


7 


7.27 


0.00 
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Table 3 
Orbital Elements 





P 


T 










f(m) 


a 1 sin i. 


a 


Star 


(day.) 


(11. ID 2. ion. nun) 






(kms 


(km s ' ) 


[M . ) 


.) 


(kms 1) 


HD 1976 


25.4176 ± 0.0004 


35783.5 ±0.1 


0.12 ± 0.03 


172 ±2 


23.6 ±0.7 


-9.7 ±0.5 


0.034 ± 0.003 


11.7 ±0.4 


4.81 


HD 14633 


15.4082 ± 0.0004 


51885.07 ±0.01 


0.700 ± 0.003 


138.7 ± 0.5 


19.0 ±0.1 


-38.17 ±0.08 


0.0040 ± 0.0001 


4.13 ± 0.04 


1.96 


HD 15137 


30.35* 


51879.7 ±0.6 


0.48 ± 0.07 


148 ± 10 


13 ± 1 


-48.4 ± 0.6 


0.004 ± 0.002 


6.7 ±0.8 


3.94 


HD 37737 


7.8400 ± 0.0002 


53690.2 ±0.1 


0.43 ± 0.02 


158 ± 5 


72 ± 2 


- 8 ± 2 


0.22 ± 0.02 


10.1 ± 0.3 


5.28 


HD 52533 


22.1861 ± 0.0002 


44635.01 ± 0.03 


0.30 ±0.01 


330.3 ±0.5 


105 ± 1 


76.7 ± 0.6 


2.30 ± 0.07 


43.8 ± 0.4 


11.7 



»Fixed. 



Table 4 

Radial Velocity Measurements of Binaries with Orbital Solutions 





HJD 


Orbital 


Vr 


(O-C) 
(km s-^) 




HJD 


Orbital 


Vr 


(O-C) 
(km s-^) 


Star 


(-2,450,000) 


Phase 


(km s-^) 


Star 


(-2,450,000) 


Phase 


(km s-^) 


rlU L\)(0 


ODO I . lol 


0.225 


— 14.2 


— 0.5 




ODyO-DZO 


0.864 


— 48.8 


— 4.4 




c; O O Q Q 
ODOO.OOO 


0.266 


—4.8 


3.0 




ODyD.o41 


0.871 


—38.3 


6.5 




rip, CO c Q c 


0.267 


— 4. .3 


3. -5 














3659.720 


0.301 


-0.3 


3.0 


HD 37737 


3658.005 


0.898 


-13.1 


5.2 




oboy. / zi 


0.301 


— 0.1 


3.2 




ODOO. woo 


0.023 


— 108.8 


—0.2 




ODOU. / o4 


0.342 


3.0 


1.7 




oDoy.yiD 


0. 142 


— 50.9 


0.2 




ODDo. ( OO 


0.459 


13.7 


4.0 




ODD.^.yoo 


0.533 


31.4 


—0.7 




ODOO. ( OO 


0.325 


— 9.4 


— 8.9 




ODDO.OiO 


0.537 


27.4 


—4.9 




ODyo.DoO 


0.636 


13.8 


6.3 




ODDo.yoO 


0.660 


30.1 


— 5.0 




obyo. ( yu 


0.642 


7.3 


0.2 




OOOO.OOO 


0.575 


27.2 


— 6.7 




ODy4. f Ui 


0.678 


6.9 


2.6 




ODO f .OoO 


0.600 


34.7 


—0.1 




ODyo.D ( 4 


0.716 


— 1.4 


— 1.8 




ODO i -ooZ 


0.702 


31.6 


—2.0 




3695. 769 


0.720 


—0.8 


—0.7 




ODOO.OoY 


0.729 


30.7 


— 1.1 




3696.505 


0.751 


—4.9 


—0.9 




oooo.yUO 


0.839 


14.5 


3.8 




Oftnfl '7QO 

3696. /o9 


0.760 


— 2.0 


3.2 




oDoy.0.^0 

m OCA 
ODyo.oD4 


0.855 

0.472 


2.6 

23.7 


— 2.0 

— 3.7 


HD 14633 


3657.799 


0.051 


-60.0 


-0.2 




3693.988 


0.488 


37.4 


8.6 




3657.803 


0.051 


-59.4 


0.3 




3694.886 


0.602 


44.8 


9.9 




3658.843 


0.118 


-50.1 


-1.1 




3695.868 


0.727 


33.4 


1.5 




3659.756 


0.178 


-40.0 


4.3 




3695.985 


0.742 


36.0 


5.5 




3660.760 


0.243 


-47.5 


-6.6 




3696.737 


0.838 


11.2 


0.3 




3663.739 


0.436 


-32.8 


2.4 




3697.012 


0.873 


-12.5 


-8.6 




3693.639 


0.377 


-36.7 


-0.1 














3693.852'' 


0.391 


-49.0 


-12.7 


HD 52533 


3658.024 


0.697 


-17.9 


-19.8 




3694.705 


0.446 


-33.4 


1.6 




3659.020 


0.742 


21.2 


11.8 




3694.863 


0.456 


-36.1 


-1.3 




3659.970 


0.785 


30.4 


8.1 




3695.679 


0.509 


-37.1 


-3.3 




3659.978 


0.786 


33.2 


10.8 




3695.875 


0.522 


-32.3 


1.2 




3662.994 


0.921 


119.4 


-0.6 




3696.624 


0.570 


-36.8 


-4.2 




3663.020 


0.923 


118.7 


-2.5 




3696.838 


0.584 


-29.6 


2.7 




3663.994 
3684.995^'^ 


0.967 
0.913 


160.6 
108.9 


-6.4 
-2.5 


HD 15137 


3657.812 


0.586 


-40.8 


0.8 






0.998 


193.8 


0.6 




3658.850 


0.620 


-45.7 


-4.5 




3687.892^''^ 


0.044 


205.8 


-2.3 




3659.806 


0.651 


-41.2 


-0.2 




3688.877"^'^ 


0.088 


169.8 


-27.0 




3659.810 


0.651 


-36.6 


4.4 




3689.006^'^ 


0.094 


181.2 


-12.8 




3660.765 


0.683 


-47.0 


-6.1 




3693.875 


0.313 


67.4 


-1.3 




3663.742 


0.781 


-40.9 


0.5 




3693.998 


0.319 


62.8 


-3.6 




3693.644 


0.766 


-39.9 


1.3 




3694.894 


0.359 


46.9 


-3.7 




3693.855 


0.773 


-47.7 


-6.4 




3695.023 


0.365 


47.8 


-0.8 




3694.713 


0.801 


-38.6 


3.2 




3695.896 


0.404 


30.4 


-5.4 




3694.876 


0.807 


-42.0 


-0.0 




3696.002 


0.409 


37.6 


3.3 




3695.683 


0.833 


-47.9 


-5.1 




3696.825 


0.446 


43.2 


18.9 




3695.887 


0.840 


-37.9 


5.2 




3697.019 


0.455 


28.1 


5.9 



^Obtained with the coude feed telescope 
^Measurement assigned zero weight 
"^Measurement assigned half weight 
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Table 5 

Radial Velocity Measurements of Other Stars 





HJD 


v; 






HJD 




a 


Star 


(-2,450,000) 


(km s"^) 


(km s"^) 


Star 


(-2,450,000) 


(km s"^) 


(km s"^) 


Feige 25 


ODo7.oo4 


34.3 


10.4 


HU o2zDb 


oDbo.yyi 


20.5 


10.9 




3bDo.oo7 


36.0 


5.6 




3684.980 


30.2 


18.4 




ODOO. / 4y 


47.3 


11.4 




ooyo.o ( 1 


39.0 


8.2 




3688.815'' 


25.3 


16.3 




3693.993 


39.2 


12.6 




3693.838 


39.4 


7.1 




3694.891 


33.7 


6.4 




3696.689* 


43.3 


12.9 




3695.019 
3695.893 


33.5 
32.7 


5.8 
9.1 


HD 30614 


3657.989 


8.9 


12.5 




3695.999 


29.5 


13.9 




3658.906 


9.0 


6.3 




3696.822 


31.7 


13.4 




3659.010 


3.2 


6.8 




3697.016 


38.8 


9.4 




3659.899 


22.7 


10.2 












3659.901 


10.5 


8.6 


HD 60848 


3658.032 


13.5 


14.1 




3662.980 


22.5 


7.9 




3659.026 


20.7 


14.9 




3663.976 


13.3 


16.2 




3659.987 


18.9 


17.8 




3686.859'^ 


2.2 


11.4 




3659.995 


26.6 


15.0 




3693.858 


2.9 


7.4 




3663.002 


19.7 


13.6 




3693.983 


1.6 


9.7 




3663.999 


21.4 


14.3 




3694.879 


22.7 


9.1 




3685.004" 


33.8 


15.4 




3695.000 


13.3 


9.3 




3693.879 


10.2 


12.7 




3695.861 


18.4 


10.7 




3694.001 


26.0 


16.7 




3695.978 


15.4 


14.6 




3694.897 


26.8 


14.4 




3696.731 


9.8 


8.6 




3695.026 
3695.898 


18.5 
21.2 


13.4 
7.4 


HD 36576 


3657.996 


62.7 


28.8 




3696.005 


26.1 


9.7 




3658.978 


37.8 


6.7 




3696.827 


21.9 


14.5 




3659.911 


21.4 


7.2 




3697.021 


17.5 


10.2 




3662.983 


44.6 


7.7 












3663.985 


35.1 


8.2 


HD 188001 


3657.649 


16.9 


6.3 




3686.866" 


62.1 


23.6 




3661.609 


20.8 


7.6 




3687.913" 


74.2 


39.0 




3663.570 


21.8 


8.5 




3693.869 


23.5 


11.4 




3687.603" 


11.8 


9.4 




3693.992 


27.4 


10.1 




3693.551 


6.3 


5.3 




3694.890 


45.1 


9.5 




3694.561 


6.2 


6.5 




3695.019 


55.4 


8.7 




3695.556 


16.0 


4.3 




3695.871 


58.3 


6.7 




3696.553 


25.3 


6.1 




3695.988 


53.9 


5.4 












3696.741 


61.4 


7.2 


HD 195592 


3657.668 


-24.5 


12.2 




3697.009 


34.8 


8.9 




3657.671 
3661.615 


-23.2 
-19.3 


3.9 
7.4 


HD 52266 


3658.013 


14.9 


18.5 




3663.574 


-19.0 


7.8 




3659.013 


17.2 


6.7 




3688.620" 


-20.4 


3.2 




3659.956 


15.4 


6.7 




3693.556 


-28.7 


6.8 




3659.961 


16.5 


6.8 




3694.566 


-26.8 


5.3 




3662.998 


24.2 


7.9 




3695.561 


-32.1 


4.8 




3663.026 


12.7 


15.1 




3696.558 


-32.8 


6.2 



^Obtained using co-added spectra from 2.1 m telescope 

Obtained using co-added spectra from coude feed telescope 
^Obtained with the coude feed telescope 
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Table 6 
Stellar Physical Parameters 



star 


Spectral 


'irf-f 




V sin ? 






Old 






expected 


adcjpted 




Star 


Type 


(K) 




(km. 


i'iii - \ ) 


It 




(pc) 




11, (11,:) 


<I ipe) 


l{(*k~i'tuicc's 


HD 1976 


B5 IV 


16100 


3.8 


160 


0.13 ±0.02 


2.60 ±0.5 


184 ± 9 


190 


3.8 ± 0.2 


5.7 


240 ± 50 


9 


HD 14633 


08.5 V 


35100 


3.95 


138 


0.13 ± 0.02 


3.18 ± 0.5 


40 ± 2 


2150 


9.2 ± 0.4 


8.3 


2040 ±110 


6, 10 


HD 15137 


09.5 V 


29700 


3.50 


234 


0.43 ± 0.04 


3.18 ± 0.5 


56 ± 7 


2650 


16.1 ± 1.9 


13.2 


2420 ± 230 


6, 10 


Feige 25 


B7 


13200 


4.0 


250 


0.12 ± 0.01 


2.94 ± 0.1 


11 ± 0.1 


2800 


3.3 ± 0.03 


3.0 


2670 ± 130 


5, 9 


HD 30614 


09.5 lao 


29000 


3.0 


118 


0.33 ± 0.04 


3.21 ± 0.3 


269 ± 21 


1200 


34.7 ± 2.7 


36.8 


1240 ± 40 


3, 7, 9 


HD 36576 


B2 IVc 


26060 


3.82 


165 


0.45 ± 0.06 


3.28 ± 0.4 


211 ± 14 


310 


7.0 ± 0.5 


7.3 


320 ± 10 


1, 2, 9 


HD 37737 


09.5 III 


29800 


3.95 


182 


0.69 ± 0.04 


2.88 ± 0.2 


69 ± 4 


2300 


17.1 ± 0.9 


7.2 


1640 ± 660 


4, 6, 8 


HD 52266 


09 V 


31000 


3.75 


248 


0.31 ± 0.01 


3.32 ± 0.4 


66 ± 5 


1700 


12.0 ± 0.9 


13.4 


1790 ± 90 


6, 7 


HD 52533 


09 V 


32400 


3.95 


270 


0.21 ± 0.03 


3.17±0.5 


42 ± 1 


2000 


9.0 ± 0.2 


7.9 


1880 ± 120 


6, 7 


HD 60848 


08 V:pevar 


31600 


4.20 


163 


0.18 ± 0.02 


3.27 ±0.5 


63 ± 4 


1900 


12.8 ± 0.8 


7.2 


1480 ± 420 


6, 7 


HD 188001 


07.5 la 


34080 


3.36 


94 


0.32 ± 0.05 


3.18 ± 0.2 


99 ± 5 


3200 


34.0 ± 1.6 


20.8 


2580 ± 620 


6, 7 


HD 195592 


09.5 la 


30460 


3.19 


114 


1.17 ±0.04 


3.02 ± 0.2 


220 ± 12 


1400 


33.1 ± 1.8 


22.1 


1170 ± 230 


6, 7 



References. — (1) Abt et al. 2002; (2) Chauville et al. 2001; (3) Crowther et al. 2006; (4) Leitherer 1988; (5) Martin 2006; (6) Martins et al. 2005; (7) Mason et al. 
1998; (8) Patriarchi et al. 2003; (9) Schaller et al. 1992; (10) van Steenberg & Shull 1988. 



Table 7 
Peculiar Space Velocities 





e 


b 


/Iq, cos S 


A*<5 


'^t.pec 




^pec 


Star 


(dog) 


(dog) 


/ 1 \ 
(mas yr j 


/ 1 \ 

(mas yr j 


/I 1 \ 

(km s j 


/I 1 \ 

(km s ) 


/I 1 \ 

(km s ) 


nu IV 1 


iio.Oo 


— iU.DO 


io. / U It u.zo 




4t.o It 1.44 


Q 9 -1- n ft 
— y.z It u.o 


±U.4t zt i.O 




1 /in 78 


1 S5 90 

— Lo.zyj 


— U.4 It U. / 


70 _|_ 1 n 


f{f^ -1- 1 9 4 
00. D It ±Z.44 


9/1 1 _|_ 1 7 


70 R -1- 1 9 PI 
r u.o It iZ.O 


nu 1010/ 




7 

— ( .Oo 


A R7 -1- n 40 
U.D ( It U.^U 


— O.Uo Hz U. i 


( .U Hz 1 ± .U 


9(^ 9 ^ A 9 


fiO 7 _|_ 1 1 Q 
OZ. / It ± ±.0 


Feige 25 


165.42 


-48.36 


4.4 ±1.1 


-0.3 +1.2 


58.7 + 21.8 


33.2 + 4.3 


67.4 + 22.2 


HD 30614 


144.07 


+14.04 


0.49 + 0.18 


7.31 + 0.45 


46.9+ 3.2 


21.9 + 2.7 


51.8+ 4.2 


HD 36576 


187.39 


-7.84 


0.15 + 1.20 


0.01 + 0.48 


6.8+ 1.9 


34.1 + 3.3 


34.7+ 3.8 


HD 37737 


173.46 


+3.24 


0.8 +0.6 


-4.2 +0.6 


25.5 + 13.9 


-12.9 + 3.1 


28.6 + 14.3 


HD 52266 


219.13 


-0.68 


-0.89 + 0.78 


-0.70 + 0.61 


18.9+ 8.5 


-4.5 + 3.1 


19.4+ 9.0 


HD 52533 


216.85 


+0.80 


-0.67 + 2.54 


-0.04+1.81 


13.2 + 27.8 


45.1 + 2.0 


47.0 + 27.9 


HD 60848 


202.51 


+17.52 


-3.43 + 0.75 


-1.41 + 0.50 


12.7+ 7.4 


2.2 + 5.7 


12.9+ 9.4 


HD 188001 


56.48 


-4.33 


-0.17 + 0.48 


-10.45 + 0.34 


81.5 + 22.0 


-2.3 + 4.3 


81.5 + 22.4 


HD l!)r)r,<)2 


82.. SO 


+2.!)() 


-2.4G ± 0.48 


1.62 ± ().r)l 


37.2+ 7.9 


-20.1 ± 2.1 


42.3+ 8.2 



